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1 Introduction 


In studies of high-energy neutrino propagation through matter two completely 
different approaches are in use: the quantum-mechanical MSW formalism [1] 
and the classical transport theory [2]. The hrst approach accounts for the 
neutrino mixing and coherent interference of propagating neutrinos (that is 
refraction) neglecting the loss of coherence due to inelastic scattering. The 
second one provides the means for estimating neutrino absorption, produc¬ 
tion and energy loss in inelastic charged and neutral current interactions, but 
is inapplicable for the inclusion of quantum effects of neutrino mixing and 
refraction. In other words, the necessary conditions for the validity of these 
approaches are mutually exclusive. At the same time there is a multitude of 
physically feasible and interesting contexts where the effects of neutrino mix¬ 
ing, refraction and inelastic collisions may be equally important or comparable. 
Well-known examples are the propagation of high-energy neutrinos produced 
in the sun from cosmic ray interactions and from annihilations of heavy dark 
matter particles captured by the sun (see, e.g., ref. [3] and references therein). 
The so-called hidden sources (the astrophysical neutrino sources opaque for 
electromagnetic radiation) [4] hold promise as the splendid “testing laborato¬ 
ries” for the effects under discussion because some types of these objects (e.g., 
a Thorne-Zytkow star [5]) have relevant dimensions and density prohles. 

The necessity of a unihed description of mixed neutrinos in hot media (like the 
CP symmetric plasma of the early Universe, a supernova core or a protoneu¬ 
tron star) has long been realized [6] and the rigorous quantum kinetic theory 
has been developed [7]. Q However, the direct application of that theory to 
the problem of high-energy neutrino transport in a normal cold matter is very 
involved and the analytic treatment of the quantum kinetic equations [7] is 
hardly possible, except for very particular cases. 

Since an analytic analysis is always useful to gain an insight into a sophisti¬ 
cated problem, it seems instructive to apply the conventional MSW formalism 
straightforwardly generalized by incorporating the imaginary parts into the 
neutrino indices of refraction. This will allow an understanding of the impact 
of the neutrino absorption in conjunction with mixing and refraction, still 
disregarding the neutrino energy loss from the neutral current interactions 
[2,10], charged current induced chains like —>• rX, r —>• [11] and 

neutrino flavor changing reactions We” —>■ etc. [12]. The 

present paper is concerned with some surprising phenomena which follow from 
this simple approach. For further simplihcation we limit ourselves to studying 


^ For further examples of extensive study of the interplay between neutrino oscilla¬ 
tions and neutrino scattering see ref. [8]. Another closely related topic is a combined 
treatment of neutrino mixing and decay [9]. 
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two-flavor case. Some general results relevant to the subject were discussed 
previously in ref. [13]. 


2 Master equation 


The generalized MSW equation for the time-evolution operator 

/c 


m = 


Saait) Sapit) 
\Spa{t) Spp{t)^ 


of two mixed stable neutrino flavors and up propagating through a matter 
with absorption can be written as 


z—S{t) = [VHoV^ + W(f)] S(f), (S(0) = 1). 


( 1 ) 


Here 


V = 


cos 9 sin 9 
- sin 9 cos 9, 


is the vacuum mixing matrix that relates flavor to mass eigenstates, 9 is the 
mixing angle (0 < 6* < 7r/2), 


Hn 


El 0 
0 Eo 


, W(f) = -p, 


- 1 0 
^ 0 np{t) - 


( 2 ) 


= \Jpi + — Pu + T^j/‘2Pu and rrii are the energies and masses of the neu¬ 

trino mass eigenstates, respectively,!^ is the complex index of refraction, 

, , 27rNr,p(t) . 5 -^ 

= 1 +- 

Pi k 


where Nq = 6.022 x lO^^cm”^, is the amplitude for the zero-angle 

scattering from particle k {k = e,p,n,...), pit) is the density of the matter 
(in g/cm^) and Yk{t) is the number of particles k per amu in the point t of the 
medium. From the optical theorem (see, e.g., ref. [14]) we have Im [fu^k{0)] = 
(py/47r)(T*°)i, {pu), where {p^) is the total cross section for z/„/c scattering. 

We assume as usual that pf, E'^ ^ max (m^). 
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This implies that 


pjm [n„(t)] = Y, >fc(t)cT*°\ {pu) = (ty 

where Aa{t) is the mean free path of neutrino Va in the point t of the medium. 

It is useful to transform eq. (1) into the one with a traceless Hamiltonian. For 
this purpose we dehne the matrix 

S(t) = exp j\v [Ho + W(F)] dt j S(t). (4) 

After substituting eq. (4) into eq. (1), we have 

iAs(() = H(«)S(«), S(0) = 1. (5) 

Here 


H(«) 



( 6 ) 


2 2 
fri ^ _ ‘ryi ^ 

Ac = A cos 26, As = A sin 26, A = — - -i-, 

(lit) = quit) + iqiit) = -py [n^H) - nc(t)]. 

The Hamiltonian for antineutrinos is of the same form as eq. (6) but one must 
keep in mind that Re[/77„fc(0)] = -Re [/^^^(O)], while and ay'^kipu) 

are different in magnitude even at very high energies. 

The neutrino oscillation probabilities are just the squared absolute values of 
the elements of the evolution operator S(t), 

pyM ^ i^a'it)] = = \s^>o.{t)f. (7) 

Taking into account eqs. (2), (3), (4) and (7) yields 

Paa'it) = A{t) Sa'ait) , (8) 

where 

,,. \ y dt' ] 1 1 r 1 1 ■ 

(t) - exp I ^ - 2 \ ■ 
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Owing to the complex potential q, the Hamiltonian (6) is non-Hermitian and 
the new evolution operator (4) is nonunitary. As a result, there are no con¬ 
ventional relations between the four probabilities given by eq. (8). Since 


Qiit) 




1 

Aa{t) 


the matrix H(t) becomes Hermitian when = Ap. If this is the case at any t, 
eq. (5) reduces to the standard MSW equation and inelastic scattering in the 
medium results in the common exponential attenuation of the survival and 
transition probabilities. From here, we shall consider the more general and 
more interesting case, when 

The extreme example of such a case is provided by mixing between the ordi¬ 
nary (a = e,/i or r) and sterile (/? = s) neutrino flavors. Since Ag = 00 , we 
have A = 2Aa and qi = —1/4Aq. So qi is nonzero at any energy. Even without 
solving the evolution equation, one can expect the penetrability of active neu¬ 
trinos to be essentially modihed in this case because, roughly speaking, they 
spend a certain part of life in the sterile state. 


Other important examples are the z/g — i't and — Ur mixings. The total cross 
sections for e or /i production are well above the one for r production within 
a wide energy range [15]. This is because of large value of r lepton mass rrir 
which leads to high thresholds and sharp shrinkage of the phase spaces for the 
charged current Ut-N reactions as well as to the kinematic correction factors 
(oc m^) to the nucleon structure functions; the corresponding structures are 
negligible in the case of electron production and small for muon production. 
The neutral current contributions are, of course, cancelled out from qj. Thus, 
in the context of the master equation (5), Ur can be treated as (almost) sterile 
within the energy range for which ^ 

A similar situation, while in quite a different and narrow energy range, holds 
in the case of mixing of z/g with some other flavor. This is a particular case for 
a normal C asymmetric medium, because of the W boson resonance formed 
in the neighborhood of E))®® ^ 6.33 PeV through the reactions 


z/gC W hadrons and z/gC W I'ii {i = e,iJ,,T). 

Just at the resonance peak, Ri 250 (see ref. [12] for further details 
and references). 

For E,y min (where mk,w,z denote the masses of the scatterers 

k and gauge bosons) and for an electroneutral nonpolarized cold medium, the 
real part of the potential q is energy independent. In the leading orders of the 
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standard electroweak theory it is [16] 


IVoYpp 

for a 


for a 

ivo (i; - V) p 

for a 

VYnP 

for a 


e and /? = /i or r, 
fi and P = T, 
e and P = s, 
fi or T and P = s, 


( 9 ) 


where 


Vo = V2 GfNo ~ 7.63 x 10"^^ eV, (Lq = ^ ~ 1.62 x 10^ km) 

V Vo / 


ttr = 


Sar^- [ln(l/r,-) — 1] 
dTT sin^ 9w 


~ 2.44 X 10 


-5 


br = 


ln(l/r^) - 2/3 
ln(l/rT-) — 1 


~ 1.05, 


Gp is the Fermi coupling constant, a is the hne-structure constant, 6^ is the 
weak-mixing angle and = {mr/mwP- Thus, for an isoscalar medium \qji\ 
is of the same order of magnitude for any pair of flavors but — z/,.; the 
ratio is about —5 x 10“^. For certain regions of a neutron- 

rich medium the value of may become vanishingly small. In this case, 

the one-loop radiative corrections must be taken into account. For very high 
energies, all the expressions (9) have to be corrected for the gauge boson 
propagators and strong-interaction effects. 


One can expect \qji\ to be either an energy-independent or decreasing func¬ 
tion for any pair of mixed neutrino flavors. On the other hand, as noted above, 
there are several cases of much current interest when \qj\ either increases with 
energy without bound (mixing between active and sterile neutrino states) or 
has a broad or sharp maximum (as for — or Fg—mixings, respectively). 
Numerical estimations, performed using the results of refs. [12] and [15], sug¬ 
gest that for every of these cases there is an energy range in which qp and qi 
are comparable in magnitude. Considering that both qp and qi are propor¬ 
tional to the density and besides are dependent upon the composition of the 
medium there may exist even more specihc situations, when \qp\ r\j I I r\j |A| 
or \qp\ |Ac| and |g/| ~ If this is the case, the refraction, absorption 
and mixing become interestingly superimposed. 


3 Eigenproblem and mixing matrix in matter 


The matrix (6) has two complex instantaneous eigenvalues, e(t) and —e(t), 
with e = ep + iej satisfying the characteristic equation 

= (q- q+) (q - q -), (lo) 
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where q± = AciiA^ = Ae^^*®. The solution to eq. (10) for er and Ej may be 
written as 

4 = ^ (^0 - ?/) + + (11a) 

Si = ^ (provided 7 ^ A^), ( 11 b) 


with 


£0 = i/A 2 - 2Acg/J + q\ > |A^|. 

The sign of Er is a matter of convention. We dehne sign {er) = sign(A) = C,. 
At that choice e = A for vacuum (g = 0) and e = if Qi = 0. 

Let us consider the behavior of Er and Ej in the vicinity of the MSW resonance, 
defined by the condition qR = qRit^,) = Ac (we suppose for simplicity that the 
matter density and composition are continuous functions of t). The accurate 
passage to the limit in eqs. ( 11 ) gives 

^ ^s\/max(l-A 2 /A 2 , 0 ), 
lim £/ = ±CA/^max(l - A2/Af,0), 

where A/ is taken to be the value of qi in the MSW point (that is A/ = qi(t-^)). 
Therefore the resonance value of IsrI (which is inversely proportional to the 
neutrino oscillation length in matter) is always smaller than the conventional 
MSW value [A^l and vanishes if A ? < a; (e, remains hnite in this case). In 
neutrino transition through the region of resonance density p = p{t^), Ei un¬ 
dergoes discontinuous jump while Er remains continuous. The corresponding 
cuts in the q plane are placed outside the circle |g| < |A| as is shown in £g. 1 . 
If Aj > A^, the imaginary part of e vanishes while the real part is kept finite. 


A distinctive feature of eq. (10) is the existence of two mutually conjugate 
“super-resonance” points q± in which e vanishes giving rise to the total degen¬ 
eracy of the levels of the system under consideration. Certainly, the behavior of 
the system in the vicinity of these points must be dramatically different from 
the conventional pattern. As noted in sec. 2 , the “super-resonance” conditions 
are physically realizable for various meaningful mixing scenarios. 

In order to simplify the solution to the eigenstate problem we will assume 
that the phase trajectory q = q(t) does not cross the points q± at any t. 
Within the framework of the non-Hermitian quantum dynamics (see, e.g., ref. 
[17] and references therein) one has to consider the two pairs of instantaneous 
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Fig. 1. Zeros and cuts of e in the q plane for Ac > 0. 
eigenvectors \^±) and \^±) which obey the relations 

= ±e\^±) and = ±e*|#±). (12) 

For q ^ q± these pairs form a complete biorthogonal and biorthonormal set, 

(1F±|!F±) = i, = \^+){W+\ + \^_){W_\ = 1. 

Therefore, the eigenvectors are dehned up to a gauge transformation 
\'P±) ^ |^±) e"*-^±|#±), 

with arbitrary complex functions f±{t) such that Im (/±) vanish as g = 0.[^ 
Thus it is sufficient to hnd any particular solution of eqs. (12). Taking into 
account that = H*, we may set |lF±) = \\P^) and hence the eigenvectors 
can be found from the identity 


H = £|iF+)(!F;| 


T 

Setting |!^±) = (n±, ±ni|:) we arrive at the equations 


v± = 


e±{q- Ac) 
2e 


A, 


a particular solution of which can be written as 

Av-^)I2^ n. = C 


v+ = 


e + q- A, 


2e 


e - q + Ac 


2e 




( 13 ) 


( 14 ) 


® For our aims, the class of the gauge functions may be restricted without loss of 
generality by the condition /±|^=o ~ 



















where 




arg(e + q — Ac) = — arg(e — q + Ac) = arctan 




arg(e) = arctan 



1 



•> 


and we have hxed the remaining gauge ambiguity by a comparison with the 
vacuum case. 


It may be sometimes useful to dehne the complex mixing angle in matter O 
by the relations sin 0 = and cos 6> = or, equivalently, 

A, Ac — q 

sin26> = —, cos26> =-. 

e e 

The real and imaginary parts of 0 are found to be 


Re(6>) = - arctan 
Im(6)) 




{qi — As) er — {qn — Ac) Si 

(qR — Ac) Er + {qi — As) El 

_L 


{qR - Ac)^ + {qi - As)^J 


Having regard to the above-mentioned prescription for the sign of Er, one can 
verify that 6> = 6^ if g = 0 (vacuum case) and 6> = 0 if A^ = 0 (no mixing or 
rrii = m^)- It is also clear that 0 becomes the standard MSW mixing angle 
with Im(6>) = 0 when qi = 0 {A^ = H/j). 


In order to build up the solution to eq. (5) for the nondegenerate case one has 
to diagonalize the Hamiltonian H. Generally a non-Hermitian matrix cannot 
be diagonalized by a single unitary transformation. But in our simple case 
this can be done by a complex orthogonal matrix (extended mixing matrix in 
matter) 


U/ = U exp(ff). 


where f = diag (/_, /+), f± are the arbitrary complex phases mentioned above 
and 


u = (|i^_>,|s^+)) 




COS0 sin6) 
sin 0 cos 0 


The matrix U by its construction has the following properties: 

U^HU = diag {-E, e) , U^U = 1, = V. 


(15) 
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From eq. (10) it follows that de/dq = {q — Ac)/e and thus eqs. (13) yield 


9v± ^ 

dq 2^2 

We therefore have 


where 


= -Q 


0 -i 
i 0 


= —hlcr 


2e^ A dt \q — q- 


According to eqs. (15) and (16), 

UjHU^ = diag (-£, e), UjU^ = 1, 
iV]V, = -fie-^'o-^e" - f. 


( 16 ) 


(17a) 

(17b) 


4 Adiabatic solution 


With the mixing matrix U/ in hand we can write the formal solution to eq. 
(5) in the most general form: 

S(() = U,(() exp [-»(«)] X,(t)Uj(0). (18) 

Here Xf(t) is an unknown matrix, $(t) = diag (—^(t), ^(t)) and ^(t) = 
is the complex dynamical phase, dehned by 

^nit) = [ eR{t')dt', ^i{t) = [ ej{t')dt'. 

JO JO 

Substituting eq. (18) into eq. (5) and using eqs. (17), we hnd 


*X/(t) 


fi(t)e-*^WF(t)e*^W + f(^)l Xf{t), 


X/(0) = 1, 


where 


F(t) = = 



It can be proved now that the right side of eq. (18) is gauge-invariant i.e. it 
does not depend on the unphysical complex phases f±(t). This crucial fact is 


10 



closely related to the absence of the Abelian topological phases in the system 
under consideration.!^ Finally, we can put /± = 0 in eq. (18) and the result is 


S(t) = U(t) exp [-2$(t)] X(t)U^(0), 
zX(t) = fl(t)F(t)X(t), X(0) = 1. 


(19a) 

(19b) 


These equations, being equivalent to eq. (5), have nevertheless a restricted 
range of practical usage on account of poles and cuts as well as decaying and 
increasing exponents in the “Hamiltonian” HF. However, the adiabatic theo¬ 
rem of Hermitian quantum mechanics (see, e.g., ref. [19]) can straightforwardly 
be extended to eq. (5) under the following requirements: 

(a) the potential g is a sufficiently smooth and slow function of f; 

(b) the imaginary part of the dynamical phase is a bounded function i.e. 
limt^oo |^/(^)| is hnite; 

(c) the phase trajectory q = q{t) is placed far from the singularities for any t. 

The hrst requirement breaks down for a condensed medium with a sharp 
boundary or layered structure (like the Earth). If however the requirement (a) 
is valid inside each layer (tj, tj+i), the problem reduces to eqs. (19) by applying 
the rule 


S(t) = S(t,0) = S(t,t„)...S(t2,H)S(H,0), 


where S (tj+i, U) is the time-evolution operator for the Tth layer. 

The requirement (b) alone is not too restrictive considering that for many 
astrophysical objects (like stars, galactic nuclei, jets and so on) the density 
p exponentially disappears to the periphery and, on the other hand, ^ 0 
as p —> 0. In this instance, the function ^/(t) must be t independent for 
sufficiently large t. But, in the case of a steep density prohle, the requirements 
(a) and (b) may be inconsistent. 

The important case of violation of the requirement (c) is the subject of a 
special study which is beyond the scope of this paper. It is interesting to note 
in this connection that, in the Hermitian case, a general adiabatic theorem 
has been proved without the traditional gap condition [20]. 

Having regard to the mentioned limitations, we can expect that eqs. (19) is 

® Note that for case the construction corresponding to eq. (18) is not generally 
gauge-invariant (and thus is unphysical) even in the Hermitian case [18] due to the 
nontrivial (irremovable) topological phases. However, it becomes invariant both in 
the Hermitian [18] and non-Hermitian [13] cases after the substitution $ F, 

with r a diagonal matrix incorporating the topological phases. 
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tolerable for finding the approximate, adiabatic solution to the problem (and 
the corrections to that solution) for many cases of pragmatic interest. We 
shall presume now that the parameters of the matter vary so slowly that all 
necessary conditions do hold for 0 < t < T. Then, in the adiabatic limit, we 
can put n = 0 in eq. (19b). Therefore X = 1 and eq. (19a) yields 

Saa{t) = n+(0)n+(t)e“*'^^*^ + n_(0)n_(f)e*'^^*\ 

Sadit) = v-( 0 )n+(t)e“*‘^^*^ - n+( 0 )n_(t)e*'^(*\ 

= r;+(0)n_(t)e-*'^W - n_(0)n+(t)e*'^W, 

Sf3f3{t) = n_(0)n_(t)e“*'^(*) + n+(0)n+(f)e*'^^*\ 


Taking into account eq. ( 8 ) we obtain the survival and transition probabilities: 
Paa{t) = A{t) I ^ sin^ [^R{t) - V5+(t)]|, 

Pai3{t) = A(t) |[/■(^)e'^^^*) - /^(t) sin^ [^^(t) - 93-(t)]|, 

P^a{t) = A(t) |[/+(t)e'^^^*) - + /^(t) sin^ [^^(t) + 99-(t)]|, 

Pppit) = A{t) I [/r(t)e'^^(*) + /+(t)e“‘^^Wj^ _ gin 2 ^ (^+(t)]|, 

( 20 ) 


where we have denoted for compactness 

P,'it) = \v,{0)v,>{t)\, ?,<t' = ±. 




e{0)e{tW 


Z>±it) = 2 [<^( 0 ) ± 'd(f)] 


In the event that the conditions 

< 4eo(t) ( 21 a' 

and 


1 1 

Affit) A^it) 


t min [ 2 lQ,(t), 2 l^(t)] 


( 21 b) 


are satished for any f G [0, T], the formulas (20) reduce to the standard MSW 
adiabatic solution 


Paait) = Ppnit) = ^ [1 + Jit)] - Ilit) sin^ [^o(t)], 

X ' 

Papit) = Ppait) = - [1 -Jit)] + Ilit) sin^ [^o(t)], 

z ^ 


(MSW) (22) 
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where 


J{t) = 

m 


- Ac [gij(O) + quit)] + gR(0)g/j(t) 


A? 


£o(0)£o(^)’ 


£o(0)£o(^) 

= [ eo{t') dt'. 
Jo 


Needless to say either of the conditions (21) or both may be violated for 
sufficiently high neutrino energies and/or for thick media, resulting in radical 
differences between the two solutions. These differences are of obvious interest 
to high-energy neutrino astrophysics. 

It is perhaps even more instructive to examine the distinctions between the 
general adiabatic solution (20) and its “classical limit”Q 


Paa{t) = exp 


= exp 


d dt' 


Jo Aait') 
A dt' 


Jo A^(f). 


' 

Pap{t) = 0, 
Ppa{t) = 0, 


(A, = 0) 


(23) 


which takes place either in the absence of mixing or for = m^. 

In the next section we consider some features of the solution (20) for media 
with constant density and composition (g = 0). For this simple case, the adi¬ 
abatic approximation becomes exact and thus free from the above-mentioned 
conceptual difficulties. 


^ Considering that ff oc A^, the solution (23) is exact (for A^ = 0) and can be 
derived directly from eq. (5). To make certain that the adiabatic solution has correct 
classical limit, the following relations are useful: 


lim e(t) = CCr b(i) - Ac] and lim \v±{t)f 
—^0 —^0 


2 (CCh ± 1), 


where Cr = sign [qR{t) - Ac]. 
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5 Matter of constant density and composition 


For definiteness sake we assume (and thus qj < 0) from here. The 

opposite case can be considered in a similar way. Let us denote 


II = 


1 A , el + qj-Al\ f fel + qY 


1 + 


L = 


4 + ^/ 


± 


TT 

ki?l 


and ^ = 


2 

Qr - A, 


^‘r + 


^R 


As is easy to see, 


4 = 


t if sign {qR - A^) = +C, 
Jzp if sign {qn - AJ = 


4 = 4 = 


^-^-=2 = 


A. 


2 e 


and sign(<p) = —C,. 


By applying these identities and formulas from sec. 4, the neutrino oscillation 
probabilities can be written as 

P„.(*) = + l¥>l) . (24a) 

Paa(«) = (/_e-‘P4+ ^ _ /2g-i/a j,j^2 _ (24b) 

P„,2(«) = P,2„(i) = + r-e-'/'' sin^ . (24c) 


5.1 Case \q\ > jA* 


Let us examine the case when d+ and d_ are vastly different in magnitude. 
This will be true when Ap 3> A^ and the factor ^ is not too small. The 
second condition holds if qR is away from the MSW resonance value A^ and 
the following dimensionless parameter 


A, 

X = — ^ 0.033 X sin 26* 

m 


Am2 \ ^100GeV\ fVo 


10-3 eV 2 / V E, 


\q.\. 


is sufficiently small. In fact we assume |x| < 1 and impose no specihc restric¬ 
tion for the ratio qR/qi- The assumption spans several possibilities: (i) small 
Aw?, (ii) small mixing angle, (hi) high energy, (iv) high matter density. The 
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last two possibilities are of special interest because the inequality |x| < 1 may 
be fulfilled for a wide range of the mixing parameters Am^ and 6 by changing 
El, and/or p. In other words, this condition is by no means artificial or too 
restrictive. 

After simple while a bit tedious calculations we obtain 

e = 1 - + C> (x^) , = K^ + 0 (x=^) , 

I+ = l + 0 (x^) , I_ = -^x^ + 0 (x^) . 

Since has been assumed to be small compared to Ap, we have 

/ x^ 

A ^ 2Aa, A+ ^ ^ T 

In general the oscillation length in matter L has no Taylor expansion over 
the parameter x and may vary within a broad range depending on other 
circumstances. For example, L ^ Tr/lgl if |g| 3> |A| and L ^ j— if 
IqiI ^0 > |As|. Due to the wide spread among the length/time scales A-j-, 
A and L as well as among the amplitudes 4 and /, the regimes of neutrino 
oscillations are quite diverse for different ranges of variable t. 

Some of the essential features are illustrated in hgs. 2 and 3 by the example 
of — Vg oscillations in an isoscalar medium {Yp = Yn = 0.5). For such a 
medium the picture of the — Vg — ^s) oscillations is the same (nearly the 
same). In all examples, the mixing parameters are taken to be Am^ = 10“^ 
eV^ and 6 = tt/4. In order to estimate the charged and neutral current 
total cross sections we apply the results of ref. [12] based on the CTEQ4-DIS 
parton distributions [21]. The z/^e contribution is neglected. Figure 2 depicts 
the survival and transition probabilities for E,y = 1 TeV and p = 0.2g/cm^ 
(x ^ 0.131 in this case). The classical limit of (that is the penetration 
coefficient for the unmixed muon neutrino flux) is also shown for comparison. 
In Eg. 3 we present the same probabilities but for = 100 TeV and (in order 
to demonstrate the density impact) for the two values of p, 10“^ g/cm^ (upper 
panel) and 3 x 10“^g/cm^ (lower panel). In these cases x 0.261 and 0.870, 
respectively. With reference to hgs. 2 and 3, one can see a regular gradation 
from slow (at t < A^) to very fast (at t > A^) neutrino oscillations followed 
by the asymptotic nonoscillatory behavior 


^ ^ Aa, A_ ^ Aa > Aa. 
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Probabilities 



t/A 


Fig. 2. Survival and transition probabilities for the Ug oscillations at Ei, = 1 

TeV and p = 0.2 g/cm^. The penetration coefficient exp(—for the case of 
unmixed muon neutrinos is shown by dashed curve. 
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Fig. 3. The same as in fig. 2 but for Ei, = 100 TeV and for two values of the matter 
density p: 10“^ g/cm^ (upper panel) and 3 x lO”'^ g/cw? (lower panel). 
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Pss{t) ~ e P/,^(t) = P^^(t) ~ 


The latter regime is the most remarkable. Since /!_ 3> /l^, a small part of 
the flux of muon neutrinos escapes through the medium of such a huge thick¬ 
ness that would be absolutely opaque for them in the absence of mixing. At 
t > (10 — 20)A^ the deviation of the survival probability P^^ from the clas¬ 
sical expectation grows exponentially. The asymptotic value of the transition 
probability P^^ is a factor of 4/x^ in excess of P^^. This may be potentially 
beneficial for the observational neutrino astrophysics, considering that after 
leaving the medium surrounding an astrophysical neutrino source, the sterile 
neutrinos may travel a very long distance in vacuum giving rise (through the 
vacuum oscillations) to a detectable amount of muon neutrinos. 


5.2 Degenerate case 


Our consideration must be completed for the case of degeneracy. Due to the 
condition qj < 0, the density and composition of the “degenerate environment” 
are hue-tuned in such a way that q = q-c_ = Ac — i |As|. The corresponding 
formulas can be derived by a formal passage to the limit in eqs. (24). A more 
simple way is however in coming back to the master equation (5). Indeed, in 
the limit of g = g_(^, the Hamiltonian (6) reduces to 


H 


I A. 




|A,|h^. 


Considering that = 0, we promptly arrive at the solution of eq. (5): 


S(t) = 1 — it |As| 


and thus, taking into account eq. (8), we have 

P„c.(t) = (l-|A,|t)'e-*/^, 

P^^(^) = (l + |A,|^)^e-*/^ 

Pag{t) = Pi3a{t) = {A,tfe-^/^. 

Since 1 /A^ = I/Aq —4 |As|, the necessary condition for the total degeneration 
is 4Aa |As| < 1 and thus 1/A = 1 /Aq, — 2 |As| > 2 |As|. The equality only 
occurs when is sterile. Figure 4 illustrates just this particular case. For 
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Fig. 4. Survival and transition probabilities for Va ^ oscillations in the case of 
degeneracy {q = q-()- The standard MSW probabilities (dotted and dash-dotted 
curves) together with the penetration coefficient for unmixed Ua (dashed curve) are 
also shown. 

comparison, the standard MSW solution (see eq. (22)) 


Paa{t) = Pss{t) = ^ [1 + cos {2A,t)] , 
1 

Pas{t) = Psa{t) = - [1 - COS (2A^t)] , 

z 


(MSW) 


is also shown as well as the classical penetration coefficient exp {—t/Aa) (with 
l/yl„ numerically equal to 4 |As|) relevant to the transport of unmixed active 
neutrinos through the same environment. 


6 Conclusions 


We have considered, on the basis of the MSW evolution equation with complex 
indices of refraction, the conjoint effects of neutrino mixing, refraction and ab¬ 
sorption on high-energy neutrino propagation through matter. The adiabatic 
solution with correct asymptotics in the standard MSW and classical limits 
has been derived. In the general case the adiabatic behavior is very different 


19 




















from the conventional limiting cases. 

A noteworthy example is given by the active-to-sterile nentrino mixing. It has 
been demonstrated that, nnder proper conditions, the survival probability of 
active neutrinos propagating through a very thick medium of constant density 
may become many orders of magnitude larger than it would be in the absence 
of mixing. The quantitative characteristics of this phenomenon are highly 
responsive to changes in density and composition of the medium as well as to 
neutrino energy and mixing parameters. Considering a great variety of latent 
astrophysical sources of high-energy neutrinos, the effect may open a new 
window for observational neutrino astrophysics. 
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